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Abstract 

A theory has been developed to explain the anomalous behavior of the magnetic sus- 
ceptibility of a normal metal-superconductor (NS) structure in weak magnetic fields at 
millikelvin temperatures. The effect was discovered experimentally by A.C. Mota et al 
jlUj . In cylindrical superconducting samples covered with a thin normal pure metal layer, 
the susceptibility exhibited a reentrant effect: it started to increase unexpectedly when 
the temperature lowered below 100 mK. The effect was observed in mesoscopic NS struc- 
tures when the N and S metals were in good electric contact. The theory proposed is 
essentially based on the properties of the Andreev levels in the normal metal. When the 
magnetic field (or temperature) changes, each of the Andreev levels coincides from time to 
time with the chemical potential of the metal. As a result, the state of the NS structure 
experiences strong degeneracy, and the quasiparticle density of states exhibits resonance 
spikes. This generates a large paramagnetic contribution to the susceptibility, which adds 
up to the diamagnetic contribution thus leading to the reentrant effect. The explanation 
proposed was obtained within the model of free electrons. The theory provides a good 
description for experimental results |10j . 



1 Introduction 



Mesoscopic systems can exhibit surprising properties at com- 

paratively low temperatures. For pure normal metals there is a length 
scale = hVp/kBT {Vp is the Fermi velocity, T is the temperature, 
ks is the Boltzmann constant) which has the meaning of a coher- 
ence length in a system with a disturbed long-range order. When this 
length is comparable with the characteristic dimensions of the system, 
the interference effects can come into play. Theoretically this was first 
demonstrated by Kulik ||4j for a thin-wall normal pure-metal cylinder 
in the vector potential field. It appears that the magnetic moment of 
such a system is an oscillating function of the magnetic flux through 
the cross-section of the cylinder, the oscillation period being equal to 
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the flux quantum of the normal metal he/ e. The effect is generated by 
quantization of the electron motion and due to the sensitivity of the 
states of the system to the vector potential fleld (Aharonov-Bohm ef- 
fect [5]). Bogachek and this author showed the existence of oscillating 
component with the period hc/e in the magnetic moment of a singly 
connected normal cylinder in a weak magnetic fleld. Oscillations with 
this period are produced by the magnetic surface levels of the cylindri- 
cal sample in a weak magnetic fleld |E| • The effect of flux quantization 
in a normal singly connected cylindrical conductor was flrst detected 
experimentally in 1976 by Brandt et al. when they were investigating 
the longitudinal magnetoresistance in pure Bi single crystals |7j,|Bl- 
This was actually the flrst observation of the interference effect of flux 
quantization in nonsuperconducting condensed matter. 

Recent advanced technologies of preparation of pure samples have 
enabled investigation of the coherent properties of mesoscopic struc- 
tures taking proper account of the proximity effect [|^. The samples 
were superconducting Nh wires with a radius R of tens of /im coated 
with a thin layer d of high-purity Cu or Ag. The metals were in 
good contact and the electron mean free path exceeded the typical 
scale E,N- The magnetic susceptibilities of copper and silver were mea- 
sured. The breakdown fleld H\j, the supercooled fleld Hsc and the 
superheated fleld Hgh were estimated as functions of temperature and 
normal metal thickness. While continuing their experiments on these 
samples, Mota and co-workers fTDI detected a surprising behavior of 
the magnetic susceptibility of a cylindrical NS structure (A^ and S 
are for the normal metal and the superconductor, respectively) at very 
low temperatures (T < lOOmi^) in the external magnetic fleld parallel 
to the NS boundary. 

Most intriguingly, a decrease in the sample temperature below a 
certain point Tj. (at a flxed fleld) produced a reentrant effect: the de- 
creasing magnetic susceptibility of the structure unexpectedly started 
growing. A similar behavior was observed with the isothermal reen- 
trant effect in a fleld decreasing to a certain value below which 
the susceptibility started to grow sharply. It is emphasized in Ref. |IT] 
that the detected magnetic response of the NS structure is similar to 
the properties of the persistent currents in mesoscopic normal rings. 
It is assumed [9j - [12] that the reentrant effect reflects the behavior 
of the total susceptibility x of the NS structure: the paramagnetic 
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contribution is superimposed on the Meissner effect-related diamag- 
netic contribution and nearly compensates it. Anomalous behavior of 
the susceptibility has also been observed in AgTa, CuNb and AuNh 
structures jTHj . 

The possibility of the paramagnetic contribution to the susceptibil- 
ity of the NS structure needs further clarification. The NS structure 
in question is essentially a combination of two subsystems capable of 
electron exchange, which corresponds to the establishment of equilib- 
rium in a large canonical ensemble (with fixed chemical potential). 
Assume that these systems are initially isolated with a thick dielec- 
tric layer. It is known that the superconductor response to the applied 
magnetic field generates superfluid screening current near the cylin- 
der surface (Meissner effect). How does the normal mesoscopic layer 
respond to the weak magnetic field? Kulik [3] shows (see above) that 
in a weak magnetic field the magnetic moment of a thin-wall normal 
cylinder oscillates with the flux. The magnetic moment oscillations 
are equivalent to the existence of persistent current. Since the ener- 
gies of the individual states andj hence, the total energy are dependent 
on the flux, the average current is nonzero. The current state corre- 
sponds to the minimum free energy, therefore the inclusion of weak 
dissipation would not lead to the decay of the current state. When 
the N and S metals are isolated, the quantum states of the quasipar- 
ticles in the A^-metal are formed at the expense of specular reflection 
of the electrons from the dielectric boundaries. The amplitude of the 
magnetic moment oscillations in the N layer is small, which is deter- 
mined by the smallness of the parameter 1 / kpR in the problem and by 
the paramagnetic character of the persistent current (when the 

magnetic field tends to zero, the magnetic susceptibility is positive). 
Thus, in the absence of the proximity effect, the total susceptibility 
of the NS structure is only governed by the diamagnetic contribution 
of the (S-layer (the paramagnetic contribution is very small). 

When the proximity effect is present in the NS structure, we as- 
sume that the probability of the electron transit from the supercon- 
ductor to the N metal is close to unity. This significantly affects 
the properties of the NS structure. The diamagnetic response of the 
superconductor persists but new properties appear, that are brought 
about by the proximity effect. Now two kinds of electron reflection are 
observed in the normal film - a specular reflection from one boundary 
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and the Andreev reflection from other. Along with the trajectories 
closed around the cylinder circle, new trajectories appear in a weak 
fleld, which "screen" the normal metal. The new trajectories of "par- 
ticles" and "holes" conflne the quantization area of the triangle whose 
base is a part of the NS boundary between the points of at which the 
quasiparticle collides with this boundary. This area is maximum for 
the trajectories touching the superconductor. It is shown below that 
at certain values of the flux through the triangle area, the electron 
density of states experiences flux-dependent resonance spikes. Thus, 
in the presence of the proximity effect, the periodic flux-induced os- 
cillations of the thermodynamic values typical of the normal layer in 
the NS structure give way to periodic resonance spikes with a period 
equal to a superconducting flux quantum hc/2e P^J. The response of 
the normal mesoscopic layer to a weak magnetic fleld {H < lOOe) is 
paramagnetic and the susceptibility amplitude is large. The picture, 
however, changes when the quantized magnetic flux through the tri- 
angle area increases and its value divided by hc/2e starts to exceed 
the highest Andreev "subband" number. A phase transition occurs 
in a certain fleld Hr. As a result, the layer is now screened only 
by the trajectories of those quasiparticles that do not colhde with the 
superconducting boundary. Their amplitudes are rather small (see 
above) against the large diamagnetic response. We can thus conclude 
that the resonance contribution to the paramagnetic susceptibility of 
the NS structure can only appear in comparatively weak magnetic 
flelds. At this condition the reentrant effect may be generated. The 
conclusion correlates well with the experimental observations |E| - 

El- 

The origin of paramagnetic currents in structure was discussed 
in several theoretical publications. Bruder and Imry |[T7j analyze 
the paramagnetic contribution to susceptibility made by quasiclassi- 
cal ("glancing") trajectories of quasiparticles that do not collide with 
the superconducting boundary. The authors |jT7| point to a large 
paramagnetic effect within their physical model. However, their ra- 
tio between the paramagnetic and diamagnetic contributions is rather 
low and cannot account for experimental results [|9j - jHj . 

Fauchere, Belzig and Blatter [18] explain the large paramagnetic 
effect assuming a pure repulsive electron-electron interaction in no- 
ble metals. The proximity effect in the N metal induces an order 
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parameter whose phase is shifted by tt from the order parameter A 
of the superconductor. This generates the paramagnetic instabihty 
of the Andreev states, and the density of states of the NS structure 
exhibits a single peak near the zero energy. The theory in jlB] essen- 
tially rests on the assumption of the repulsive electron interaction in 
the N metal. Is the reentrant effect a result of specific properties of 
noble metals? or Does it display the behavior of any normal metal 
experiencing the proximity effect from the neighboring superconduc- 
tor? Only experiment can provide answers to these questions. We 
just note that the theories of jTTj, flB] do not account for the temper- 
ature and field dependencies of the paramagnetic susceptibility and 
the nonlinear behavior x of the NS structure. The current theories 
cannot explain the origin of the anomalously large paramagnetic reen- 
trant susceptibility in the region of very low temperatures and weak 
magnetic fields. 

It is worth mentioning the assumption made by Maki and Haas 
112] that below the transition temperature (~ lOmK) some noble 
metals (Cw, Ag, Au) can exhibit j9-ware superconducting ordering, 
which may be responsible for the reentrant effect. This theory does 
not explain the high paramagnetic reentrant effect either. 

In this paper a theory of the reentrant effect is proposed which is 
essentially based on the properties of the quantized levels of the NS 
structure. Levels with energies no more than A (2A is the gap of 
the superconductor) appear inside the normal metal bounded by the 
dielectric (vacuum) on one side and contacting the superconductor on 
the other side. The number of levels no in the well is finite. Because 
of the Aharonov-Bohn effect [5J, the spectrum of the NS structure 
is a function of the magnetic flux in a weak field. The specific fea- 
ture of the quantum levels of the structure is that in a varying field 
H (or temperature T) each level in the well periodically comes into 
coincidence with the chemical potential ( of the metal. As a result, 
the state of the system suffers strong degeneracy and the density of 
states of the NS sample experiences resonance spikes. 

It is shown that the phenomenon of resonance appears in a cer- 
tain interval of weak magnetic fields at temperatures no higher than 
a hundred of millikelvins. Resonance is realizable only in pure meso- 
scopic N layers under the condition of the Aharonov-Bohm effect. 
The resonance produces a large paramagnetic contribution to the 
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susceptibility of the NS structure. When is added to the dia- 
magnetic contribution x'^ produced by the Meissner effect, the total 
susceptibihty displays the features of the reentrant effect [20j . 

2 Spectrum of quasiparticles of the NS structure 

Consider a superconducting cylinder with the radius R which is 
covered with a thin layer d of a pure normal metal. The structure 
is placed in a weak magnetic field H{0, 0, H) oriented along the sym- 
metry axis of the structure. It is assumed that the field is weak to 
an extent that the effect of twisting of quasiparticle trajectories be- 
comes neghgible. It actually reduces to the Aharonov-Bohm effect [H], 
i.e. allows for the increment in the phase of the wave function of the 
quasiparticle moving along its trajectory in the vector potential field. 

We proceed from a simplified model of NS structure in which the 
order parameter magnitude changes stepwise at the NS boundary. It 
is also assumed that the magnetic field does not penetrate into the 
superconductor. The coherent properties observed in the pure normal 
metal can be attributed to its large "coherence" length at very low 
temperatures. 

One can easily distinguish two classes of trajectories inside the nor- 
mal metal. One of them includes the trajectories which collide in 
succession with the dielectric and NS boundaries. The quasiparticles 
moving along these trajectories have energies e < A and are localized 
inside the potential well bounded by a high dielectric barrier (~ leV) 
on one side and by the superconducting gap A on the other side. On 
its collisions, the quasiparticle is reflected specularly from the dielec- 
tric and experiences the Andreev scattering at the NS boundary jlH] . 
We introduce an angle a at which the quasiparticle hits the dielec- 
tric boundary. The angle is counted off the positive direction of the 
normal to the boundary (Fig. 1). in this case the first class contains 
the trajectories with a varying within < a < q;c (o^c is the angle 
at which the trajectory touches the NS boundary). The other class 
includes the trajectories whose spectra are formed by collisions with 
the dielectric only, i.e. the trajectories with a > a^- 

The two groups of trajectories produce significantly different spec- 
tra of quasiparticles. The distinctions are particularly obvious in the 
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Figure 1: Two classes of trajectories in the normal metal of NS structure in the magnetic field: 
a) trajectories forming the Andreev levels; b) trajectories colliding only with the dielectric 
boundary. 



presence of the magnetic field. The trajectories with a < o^c form a 
spectrum of Andreev levels which contains a supplement in the form 
of an integral of the vector potential field. The spectrum character- 
izes the magnetic fiux through the area of the triangle between the 
quasiparticle trajectory and the part of the A^^S* boundary. It is also 
determines the magnitude of the screening current produced by " par- 
ticles" and "holes" in the N layer. These states are responsible for 
the reentrant effect. The trajectories with a > ac do not collide with 
the NS boundary. The states induced by these trajectories are prac- 
tically similar to the "whispering gallery" type of states appearing in 
the cross-section of a solid normal cylinder in a weak magnetic field 
[21] • The size of the caustic of these trajectories is of the order of 
the cylinder radius, i.e. they correspond to high magnetic quantum 
numbers m. The spectrum thus formed carries no information about 
the parameters of the superconductor and it is impossible to meet the 
resonance condition in this case. These states make a paramagnetic 
contribution to the thermodynamics of the NS structure but their 
amplitude is small (~ l/kpR)- It is therefore discarded from further 
consideration. Our interest will be concentrated on the trajectories 
with a < ac- 

The spectrum of quasiparticles of the NS structure can be obtained 
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easily using the multidimensional quasiclassical method generalized 
for the case of the Andreev scattering in the system jTHj, [22] • Af- 
ter collision with the NS boundary the "particle" transforms into a 
"hole". The "hole" travels practically along the path of the "particle" 
but in the reverse direction. In the strict sense, however, the path of 
the "hole" is somewhat longer because under the condition of Andreev 
elastic scattering the momentum of the " particle" exceeds that of the 
reflected " hole" . According to the law of conservation of the angular 
momentum, the angle a' at which the "hole" comes up to the dielec- 
tric boundary and hence the distance covered by the "hole" are larger. 
Eventually, the trajectory of the quasiparticle becomes closed due to 
its displacement along the perimeter of the layer. However, as the 
quasiparticle energy decreases and approaches the value of the chem- 
ical potential, the difference a — a' starts tending to zero. Since our 
further interest is concerned with low-lying Andreev levels, we assume 
that the "hole" trajectory is strictly reversible. The distance covered 
by the "particle" ("hole") between two boundaries is jCq ~ 2d/ cos a. 

According to the multidimensional quasiclassical method fTB] , [221 ? 
there are two congruences of "particle" rays - towards the dielectric 
(/) and in the opposite direction (//). There are also two congruences 
of "hole" rays - towards the NS boundary (///) and away from it 
{IV). The covering space is constructed of four similar structures 
whose edges are joined in accordance with the law of quasiparticle 
reflection from a dielectric and a NS boundary. At the dielectric 
boundary the congruences / and // are joined. The congruences /// 
and IV are joined independently. The covering space consists of the 
outer ("particles") and inner ("holes") toroidal surfaces. Each surface 
contains only a part of the single independent integration contour. 
The path of the "particle" is 2d. The "hole" travels the same length 
whereupon the trajectory of the quasiparticle closes. The total length 
of the closed contour along the covering surface of the NS structure 
is 4d. 

It is possible to choose two independent integration contours within 
a tours that do not contract into a point. One condition of quantiza- 
tion relates the caustic radius to the magnetic quantum number m. 
We replace it with an angle of incidence of the quasiparticle on the 
dielectric boundary. The other condition of quantization introduces 
the radial quantum number n. Thus, the complete set of quantum 
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numbers describing the motion of the quasiparticle includes n, a, q, 
where q is the quasimomentum component along the symmetry axis 
of the cylinder. 

Assume that the condition d <^ R is obeyed for the NS structure. 
We can then neglect the curvature of the cylinder boundary and as- 
sume that it is flat. The condition of quasiclassical quantization can 
be written as 

Pi + —Aj ds = 27rh + 1 + — arccos e / 

Co 

(1) 

where po {pi) are the quasimomentum of the "particle" ("hole"), e is 
the "quasiparticle" energy, A is the vector potential {0,0, Hy), \jCq\ is 
the trajectory length covered by the "particle" ("hole"). The unity in 
the right-hand side of Eq. (jT]) appears when two collisions of the quasi- 
particle with the dielectric boundary are taken into account |j22] • The 
term ;^ arccos s/ A accounts for the phase delay of the wave function 
under the Andreev scattering of quasiparticles [ISj. The quasimo- 
mentum pq and pi in Eq. dH) can be expanded in the parameter £/( 
retaining the first-order terms and replacing n + 1 by n. As a result, 
Eq. dH) furnishes the sought for spectrum of the NS structure in a 
weak magnetic field {C is the quasiparticle trajectory): 

, Trhv r(q) cos a ( \ e tga^\ , , 

en{q, a- $ = n + - arccos - - . 2 

Here vc{q) = \/pp — q^ /m* , pp is the Fermi momentum, q is the 
quasiparticle momentum component along the cylinder axis, m* is the 
effective mass of the quasiparticle, $o = hc/2e is the superconducting 
flux quantum. The positive a- values refer to "particles" (n > 0), 
while the negative ones are for "holes" (n < 0). 

The last term in Eq. (I2j) has the meaning of " phase" 

d 

A{x)dx, (3) 




which is dependent on the vector potential field and varies with the 
angle a characterizing the trajectory of the quasiparticle. 

The spectrum of Eq. (|2I) is similar to Kulik's spectrum [23] for 
the current state of an SNS contact. However, Eq. (Ej) includes an 
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angle-dependent magnetic flux instead of the phase difference of the 
contacting superconductors. 

The value of the " phase" (flux) controls the diamagnetic and para- 
magnetic currents in the NS structure. To calculate it, we should 
know the distribution of the vector potential fleld inside the normal 
metal. 

The problem of the Meissner effect in superconductor-normal metal 
(proximity) sandwiches was solved by Zaikin [21] . It was shown that 
the proximity effect caused the Meissner effect bringing an inhomo- 
geneous distribution of the vector potential fleld over the layer of 

the structure: A(x) = Hx + ^j{a)x{d — |). For convenience we in- 

d 

troduce the notation a = J A{x)dx. This expression can be obtained 



from the Maxwell equation rot^ = with the boundary conditions 
A{x = 0) = and dxA{x = d) = H. The screening (diamagnetic) 
current j is a function of a, j{a) = —js'-p{a/^o), where jg is the super- 
fluid current and i^{x) is function of flux. Thus, we can write down 
the self-consistent equation for a j2SI - I2B] : 

Hd^ 47r . , , . , . 

« = — + (4) 

The diamagnetic current j^(a) was calculated in terms of the mi- 
croscopic theory as a sum of currents of quasiparticles ("particles" and 
"holes") for all quasiclassical trajectories characterized by the angles 
and ip [211, EBj (below the system of units ks = h = c = 1 is used): 



/($,r) = 

7r/2 n/2 

= -AT ^ f dO j dip- 



sin^ 6 cos p sin[2$tg6' cos p] 



Wn>0 







^^'+^' sha^ + ^cha 



2 



+ cos2($tg6' cosp) 
(5) 

where A = 2ekp/7T'^, Un = (2n+ l)7rT, 2A is the superconductor gap, 
an = 2ujnd/vF COS 9, and $ is given by Eq. (EI)- The function j^($) 
is noted for interesting features. In small magnetic flelds ($ <C 1) 
jd ^ —js^. Such low flelds can lead to the effect of extrascreening 
of the external magnetic fleld (see [21|)- When the fleld increases 
($ ~ 1), the current starts oscillating and for certain "phases" it 
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turns to zero at regular intervals "phases" With high values of the 
inequality ^ 1), the current amplitude decreases. 

3 Resonance spikes in the density of states of NS 
structure in weak magnetic fields 



In the region of weak magnetic fields, the density of states of the 
quasiparticles that are described by the spectrum of Eq. (j2I) exhibits 
sharp singularities. The spectrum of Eq. ((21) is formed by the tra- 
jectories of the quasiparticles which collide with the dielectric and 
superconducting boundaries. It encloses a certain area penetrated by 
a magnetic flux. At any instant when the magnetic flux becomes a 
multiple of the superconducting flux quantum, the density of states 
experiences resonance spikes. 

Let us consider the cross-section of a NS structure. Assume that 
the superconducting cylinder radius R and the normal layer thickness 
d have a mesoscopic scale. The density of states z/(£) can be calculated 
proceeding from the expression 

iy{e)=^ j dq5[£-£n{q,a)]. (6) 

The summation is taken over all quantum numbers n, a and 
spin a. Since we are not interested in the contribution from the states 
formed by the trajectories of the quasiparticles with a > we can 
write down 

Olc 

iy{£) = J dav{£;a), (7) 



-etc 



where z/(e; a) is the contribution to the density of states from the pre- 
assigned trajectory with a fixed a. Eq. ((21) for the low- lying Andreev 
levels (e <C A) is taken as a spectrum. After integration with respect 
to q and introduction of the notation P = 7r/l/2(im*, we can pass 
on to the dimensionless energy e = £/(3pF- For we have the 

expression 

2pF 2 sec^a + x| — eseca] , . 
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where x = 1/2 — ^iga/iT, and 6{x) is the stepwise Heaviside function. 
Eq. dHD suggests two cases depending on the parameter n + x. 

a) Non-resonance case. If n + x 7^ 0, the energy dependence under 
the radical sign in Eq. ^ can be neglected for small energies ( e ^ 0). 
Then, the nonresonance contribution to the density of states is 



(0) '^Pf 2 



da 



00 9 
sec a 



n + >cY 



(9) 



The series in Eq. (P) is calculated readily by the formula in j2Z| : 



E 



= (-1) 



n-l 



TT 



After calculation of the integral we obtain 



(0) Pf 2^0, 2 



a 



27ra /2i? 



(10) 



where y ^ ~ tgo^c- 

b) Resonance case. Now we go back to Eq. (Ej). We find v^"^^ as 



, res 




sec^a6'[|a„ — 6tga| — eseca] 
\(Xn — 6tgapA/|cin — &tgap — e^sec^a' 



'ir 



where the notations an = n+ l/2, 6=|^ are introduced. Eq. ( ITTl ) 
shows that at certain values of the flux (b), the radicand in the de- 
nominator turns to zero. 

Prior to calculation of v^^^, let us discuss the question of the con- 
tribution of different angles a to the resonance amplitude. It is rea- 
sonable to assume that because of the factor sec^a in the numenator 
of Eq. dm), the angles a ~ o^c are the main contributors to the 
integral. It is convenient to employ in the integral a new variable 
of integration x = tga. Then the neighborhood of the upper limit 
xo = tgofc is the main contributor to the integral. Introducing the 
notation a = a„ — 6x0 and the small deviation ^ = — a: <C 1, we can 
write down the equation for the resonance condition as: 



(62-e2)^2^ 2(56 + 6^x0)^ 



a 



6^(1 + ^2) =0. 



(12) 
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The point of our interest is the asymptotics z/(e) at low e ^ 0. Eq. 
( IT^ is solved to the accuracy within first-order terms of |e|: 



The expression in front of the radical in the denominator of Eq. 
( ITTl ) has the second order smallness in |e|, i.e. \a\^ > + Xq), 

which leads to its cancellation with the similar small parameter in the 
numenator. 

The remaining integral is estimated to be a constant of about unity. 
The resonance -induced spike of the density of states always appear 
when the Andreev level coincides with the Fermi energy at a certain 
flux in the layer. In the vicinity of the chemical potential there 
is a strong degeneracy of the quasiparticle states with respect to the 
quantum number q. As a result, a macroscopic number of q states 
contribute to the amplitude of the effect. Near the resonance, the 
ratio of the resonance and nonresonance amplitudes of the density of 
states is 

^res I 

It is thus shown that a change in the magnetic flux leads to reso- 
nance spikes in the density of states of the NS structure. The flux 
interval between the spikes is equal to the superconducting flux quan- 
tum $0- 

4 Calculation of susceptibility of NS contact 



To explain the reentrant effect, we need to have an expression for 
the susceptibility of the NS structure. We assume that in a weak 
magnetic field the total susceptibility of the NS sample consists of 
two contributions. Firstly, the response of the superconductor to the 
applied magnetic field generates the Meissner effect. Note that the 
diamagnetic response is observed in all fields up to the critical one. 
The amplitude of the diamagnetic current increases monotonously 
with lowering temperature. On the other hand, the presence of a 
pure normal metal in the NS structure produces a paramagnetic 
contribution. In a weak magnetic field the contribution is due to 
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the Aharonov-Bohm effect and the quantization of the quasiparti- 
cle spectrum of the mesascopic system. When the penetrabihty of 
the barrier between the metals is small, the electrons of the normal 
metal are reflected specularly from its boundaries. As compared to 
the diamagnetic contribution from the superconductor, the paramag- 
netic contribution produced by the N layer has a small amplitude and 
can therefore be neglected. Thus, the paramagnetic and diamagnetic 
contributions cannot compete in the absence of the proximity effect 
in the NS structure. However, if the penetrability of the barrier is 
close to unity, the mechanism of the Andreev reflection becomes ac- 
tive at the NS boundary The quasiparticle spectrum of the N layer 
undergoes a significant transformation and resonance spikes appear 
in the amplitude of the density of states in a certain regions of mag- 
netic fields and temperatures. Simultaneously, the distribution of the 
vector potential field in the normal layer becomes inhomogeneous. As 
shown below at certain values of the parameters of the problem, the 
paramagnetic contribution to the susceptibility of the NS structure 
can become equal to the diamagnetic contribution. This is the reason 
why the reentrant effect appears in pure mesoscopic NS structures. 

Theoretically, the resulting susceptibility including the reentrant 
effect can be represented as a sum of the paramagnetic contribution 

of the NS structure caused by the Andreev scattering and the dia- 
magnetic susceptibility of the system in which there is no proximity 
effect between the N and S metals. The temperature-induced behav- 
ior of the diamagnetic current in such a system is well known. As the 
temperature decreases, the diamagnetic current amplitude increases 
and becomes saturated at temperatures about several millikelvins. At 
high temperatures /c^T ^ hVp/d, the diamagnetic current decreases 
rapidly following the law j ~ exp{—4:7TkBTd/hVF). Note that in 
a NS structure in which the electrons are reflected specularly at both 
boundaries of the normal metal, the susceptibility is negative (i.e. dia- 
magnetic) in the whole interval of temperatures < T < Tg. However, 
we will not use this approach to estimate the resulting susceptibility. 
Below we calculate the screening current of the NS structure. It natu- 
rally allows for the paramagnetic contribution at certain values of the 
magnetic fleld and temperature. We focus our attention on calcula- 
tion of the paramagnetic contribution in structures with a pronounced 
proximity effect. This is important especially in the context of the re- 
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cent statement [2B] that no paramagnetic reentrance can occur in NS 
proximity cylinders in the absence of electron-electron interaction in 
the A/^ layer. 

a) Paramagnetic susceptibility of A^^S* contact 



The contribution of the states in Eq(l2l) to the paramagnetic sus- 
ceptibility of the normal layer in a NS contact can be calculated pro- 
ceeding from the expression for the thermodynamic potential (/c^ = 1) 

n = -r^ln[l + exp(-£„(g,a)/r)], (15) 

where the summation is taken over the spin (cr) and all the states 
related to the trajectories of the quasiparticles with a < Oc- The 
expression for susceptibility (per unit volume V of the normal metal) 
is found using the formula 

1 d^n 

After performing the summation over the spin and taking into ac- 
count two signs of the angle a and of the quasimomentum component 

we arrive at the initial expression for paramagnetic susceptibility 
{( is the chemical potential of the metal): 

oo 

d f d£exp{£/T) 
^ " 2rm*2$2 J [exp(e/r) + l]2^ 

(16) 



< 

etc Pf 



X 



J da cos a sin^ a J dq{p^p — q'^)^^^5{£ — £n{q,CY))- 



In Ref.|[2Dl we lost one of the radicals {p^p — q^y/^ in the similar 
initial expression for x- As a result, the amplitude of the paramag- 
netic contribution appeared to be underestimated. This mistake is 
corrected in this work. 

It is convenient to present the spectrum in terms of = -^^^ and 
X = ^ — as £n{q,a) = P cos a{n+ k) y/pp — q^. Now we introduce 
the dimensionless energy e = ^ = f^, 5£ = ^^^^ is the distance 
between the Andreev levels in the ^A^ structure. Since C/5£ ^ 1, 
the lower limit of the energy integral can be replaced with —oo. By 
introducing the variable x = tga and the notation a„ = n + 1/2, 
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b = b{H,T) = 2a/$o, cl = J A{x)dx, xq = tgao = \/2R/d and taking 



into account the parity of the integrand we obtain, instead of Eq. 




dee^ f dxx'^ 9[an — bx — e\/l + a:^ 




(17) 



In Eq.(ll7l) the summation is taken over the quantum numbers of 
the "particles". Here C = ^r^, r] = ^, no is the number of Andreev 
levels in the potential well and 9 is the Heaviside step function. It is 
seen in Eq.(IT7|) that for the given "subzone" n the amplitude of the 
paramagnetic susceptibility increases sharply whenever the Andreev 
level coincides with the chemical potential of the metal. The resonant 
spike of susceptibility occurs when an — bx tends to zero on a change in 
the magnetic field (or temperature). Because of the finite number of 
Andreev levels, the existence region of the isothermal reentrant effect 
is within < H < Hmax- 

Let us calculate the integral over x in Eq.([T7j). It contains a singu- 
larity under the radical R{x) = v^]4aS^~+~Ba^+~C where A = b'^ — e^, 
B = — 2a„6, C = — e^. The singularity is determined by the roots 
of the quadratic equation xi^2 = ^ri^ i ■^^\/b'^ + ~ 
troducing the notation = the expression for the roots can be 
written with a linear accuracy with respect to e as 



|e| 

The main contribution to the integral over Eq.(ll7l). is made by 
the vicinity of the point e ^ 0. If we exclude the singular points 
from the interval of integration, the indefinite integral over x can be 
calculated accurately (see the details in the Appendix). Because the 9- 
function is present under the integral, the integration intervals (0, xi) 
and (x2, Xq) make a finite contribution to the integral. On substituting 
the limits of integration, the expressions obtained have different pow- 
ers of the parameter |e|^^. We retain only the most important terms 
in order |e|~^ that determine amplitude of the effect. The discarded 
terms have higher orders of e -smallness. The intervals (0,a:i) and 
(x2, Xq) make contributions of the same order of e-magnitude. The 
region (xi, X2) does not contribute to the integral at all. 
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The estimate for the integral over x is 



3 5(l + a2)2g4- 



(19) 



On substituting Eq.( IT^ into Eq.(IT7l). the parameter drops out of 
the energy integral and we can take it quite easily. Taking into account 
the energy limits 0{an — |e|) appearing in the process of calculation 
we can obtain the expression for the paramagnetic contribution to the 
susceptibility of the NS structure, which in dimensional units has the 
form 



In Eq. ( 1201) the summation over the quantum number n is taken 
within finite limits, where no has the meaning of the maximum number 
of the Andreev levels inside the potential well of the NS structure. Its 
order of magnitude in no ~ where Se is the distance between 

the Andreev levels, Se = TrhVp/^d, and 2A is the energy gap. The flux 
T) = 2a/$o depends on both the magnetic field and temperature. 
In the pre-assigned field its value is dictated by the screening current 
of the NS structure j = —js(p{a/^o) (see Eq.dlj)). The obtained 
expression for manifests a more rapid decrease susceptibility at 
the increasing parameter b{H,T) than it was evidenced by Eq.([HI) in 



We first discuss the isothermal case of a very low temperature and 
clear up the qualitative behavior of susceptibihty in Eq.(12Dl). We shall 
proceed from the region of very strong magnetic fields (a/$o ^ 1) in 
which the second term in Eq.(|lj) is negligible. Then the dimensionless 
flux b{H, T) ^ 1 and the amplitude of the paramagnetic contribution 
in Eq. (EDI) decreases as b{H, T) raised to power 3. In comparatively 
weak magnetic fields a/$o ~ 1), the function ip{x) is actually an 
oscillating function of H and here we can expect the reentrant effect. 
Indeed as the field decreases to a certain value and the parameter 
b{H^T) /jiQ becomes ~ 1 (no is the number of the Andreev levels in 
the potential well), the amplitude of the paramagnetic susceptibility of 
the NS structure accepts for the first time an appreciable contribution 
from the highest Andreev "subband" (level). On a further decrease in 




(20) 



Ref.pg. 
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this field, the contribution from the highest "subband" persists, but 
in a certain lower field an additional contribution appears from the 
neighboring lower-lying "subband" uq — 1. Finally, in a very weak 
field all the "subbands" of the NS structure start to contribute and 
the paramagnetic susceptibility amplitude reaches its peak. However, 
at H ^ (a/$o ^ 0), the paramagnetic contribution turns to zero, 
as follows from Eq. (EDI)- The reason is that the resonance condition 
for the Andreev levels (Eq.(|2I)), cannot be realized in a zero field. 

Now we change to the case when the temperature of the NS struc- 
ture varies but the field is kept constant. We assume the field to 
be weak {H ~ 2 • lO^^Oe). The second term in Eq.dlj) for the fiux 
is very important. It is highest at millikelvin temperatures. As a 
result, the parameter b{H,T) has the lowest value. In this tempera- 
ture region the hyperbolic tangent is close to unity and the paramag- 
netic contribution is dependent only on the parameter b{H, T). Under 
this condition, all the "subbands" of the NS structure contribute to 
the amplitude of the effect. As the temperature rises, the parameter 
6(i7, T) increases smoothly. Simultaneously, the argument of the hy- 
perbolic tangent decreases. At a certain temperature, when the con- 
dition ksT > nhVp/'id is met, the contribution from the lowest "sub- 
band" starts dying down and its amplitude is decreasing linearly with 
growing T. On a further rise of the temperature, the contributions 
from the higher "subbands" of the spectrum die down in succession. 
Finally, at a very high temperature the paramagnetic contribution 
tends to zero. 

Let us estimate the amplitude of the paramagnetic contribution. 

The parameter 6(i7, T) is dependent on the value of the fiux a = 

d 

J A{x)dx, which at constant T can be found by solving the self- 



consistent equation Eq. OH). In the region of millikelvin temperatures 
and magnetic fields i7 ~ 2 • lO^^Oe the paramagnetic contribution 
has the largest amplitude. We obtain b{H, T) ~ 10^^ in this re- 
gion of T and H. The coefficient before the sum in Eq. (EDI) can 
be found by substituting ~ 8.75 • lO'^'^erg, d = 3.3 • lO^^cm 
Vp'^ ~ 1.39 • lO^cm/ sec for the characteristic parameters of the nor- 
mal Ag layer. We thus obtain IGC^d^ /STrhVp^l ~ 2.418 • 10^. The 
product of this coefficient and the parameter 6(i7, T) yields the order 
of magnitude of the paramagnetic contribution amplitude. It is seen 
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that the largest amphtude of the paramagnetic contribution exceeds 
that of the diamagnetic contribution in the vicinity of T = 0. 

b)Full magnetic susceptibility of NS structure in the presence of 
proximity effect 

Let us consider a structure in which the electrons experience the 
Andreev scattering at the NS boundary. In the presence of mag- 
netic field, the screening current is induced in the normal layer due to 
the Meissner effect. We estimate the susceptibility generated by this 
current. 

The total current J is related to the magnetic moment M as 

M = -JSo, (21) 
c 

where 5*0 c::^ ttR"^ is the cyhnder cross-section (d <^ R). Let the average 
current density be j. The total current is then J = S'j, where S = dL 
(L is the cylinder generatrix). The density of the screening current 
in NS proximity sandwiches was calculated by Zaikin [24j,[28j. We 
reproduce the formula for the current density (see Eq. ((SI), which 
is valid at arbitrary values of temperature and magnetic field. At 
T < hVp/d it is 

tt/2 7r/2 

4eklT ^ f f , . 0^ sin[2tg^cos(^$] 
TT^ J sh + cos^ tg6' COS 

(22) 

Here a„ = y^^of 6" ~ (^^ l)7rT and the phase $ follows from 
Eq.dSl)- Near T = the summation of frequencies in Eq. ( 0^ can be 
replaced with integration. For $ < 1 the response of the current is 



n/2 tt/2 

j($) ~ d6 dipsiii^ 6cos6cosipsm[2^tg6cos(p]. (23) 



If the field is small enough to meet the condition $ <C 1, Eq. 
reduces to the result that was obtained for the first time in [24j : 

At "phases" $ 1, the screening current of Eq. ( OHj ) turns to zero. 
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Figure 2: a) Dependence of screening current Eq. (5) on "phase" $ at T = 0, 092if. Current 
is in arbitrary units, b) The derivative of current with respect to "phase" (magnetic field) 
changes its sign in the region of low $-values. When the magnetic field tends to zero the full 
magnetic susceptibility of NS structure is positive. 

The current-phase dependence at T = is plotted in Fig. 2. The 
dependence is nonhnear and its amphtude has a maximum at a certain 
value of Knowing the current-phase dependence, we can determine 
the susceptibihty of the NS structure using the equation x = dM/dH. 
It is seen in Fig. 2 that the susceptibility x of the NS structure (the 
derivative of current with respect to field) changes its sign at a certain 
low value of the magnetic field Hj.. The susceptibility is " diamagnetic" 
in the region of high magnetic fields and "paramagnetic" at H < Hr- 
The " paramagnetic" portion of the curve is due to the proximity effect 
at the NS boundary and to the Andreev levels in the layer. 

Let us estimate x the linear-response regime near T = when 
this dependence is described by Eq.fE^. In such weak fields we obtain 
$ ~ — if^' where the "penetration depth" \n into the normal 
metal is dependent on temperature [26j: 



A^'(O) 



(r~o) 



(25) 




The estimate of susceptibility in the millikelvin region is x 




cm 
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R = 8.2 ■ 10-^cm, kp^ ~ 1.2 • 10^cm-\ Vp^ - 1.39 • lO^cm/sec, 
\n{^) ~ 2 • lO^^cm we obtain x = —0.06, which is close to x = 

Now we keep the magnetic field (assuming it weak) constant and 
plot the screening current versus temperature in a wide T-range. This 
dependence plotted using Eq. (EI) is shown in Fig. 3. It is seen that 
the current amplitude has a maximum at a certain T^. 

5 Discussion 

In this study we have investigated the behavior of a superconduct- 
ing cylinder covered with a thin layer of a pure normal metal. It is 
assumed that the normal metal and superconductor are in good con- 
tact. The system was placed in a magnetic field directed along the 
NS boundary. The NS structure has mesoscopic scale dimensions. 
It is assumed that the mean free path of the quasiparticles in the N 
layer exceeds the characteristic length = KVp/kBT, which has the 
meaning of the coherence length for a system with disturbed long- 
range order. The goal of this study was to interpret the experiments 
in which A.C.Mota et al. [lOj - [14J observed anomalous behavior of 
the magnetic susceptibility of a NS structure at varying temperature 
in a constant magnetic field or in a varying magnetic field at a con- 
stant temperature. This phenomenon was called a reentrant effect. 
Until recently it has not been explained adequately. 

Earlier [20j the author clarified the nature of the reentrant effect. 
As was found, the origin of the paramagnetic contribution is closely 
connected with the properties of the quantized Andreev levels that 
are dependent on the magnetic flux varying with both temperature 
and magnetic field. Typically, the levels in the NS structure time 
from time (at certain values of the field H or temperatures) coincide 
with the chemical potential of the metal. As a result, the state of 
the system is highly degenerate and the density of states of the NS 
structure experiences resonance spikes. The response of the normal 
mesoscopic layer to a weak magnetic field is paramagnetic. 

A theory of the reentrant effect has been developed in this study. 
We calculated the paramagnetic contribution separately and analyzed 
its behavior in a varying magnetic field and at varying temperature. 
In the course of this calculation we corrected the mistake made in [2II| 
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which led to underestimation of the effect amphtude. The paramag- 
netic response is determined only by the trajectories of the quasiparti- 
cles that collide with the NS boundary. It is shown that the reentrant 
effect can occur in a certain range of weak magnetic fields at tempera- 
tures no higher than lOOmK. We believe that paramagnetic reentrant 
effect is an intrinsic effect of mesoscopic NS proximity structures in 
the very low temperature limit. 

Assume that the temperature of the NS structure is about 10~^K 
and the magnetic field is increasing. As soon as the field exceeds a cer- 
tain value Hr, the isothermal reentrant effect must vanish. In strong 
fields the Andreev levels cease to make a resonance contribution to 
the paramagnetic susceptibility. Now the paramagnetic contribution 
is made by the states formed by the trajectories of the quasiparticles 
that collide only with the dielectric boundary. However, their contri- 
bution to the resulting susceptibility of the structure is small because 
of the smallness of the quasiclassical parameter of the problem 1 / kpR. 
Under this condition the susceptibility exhibits diamagnetic behavior 
in all strong fields up to the critical one. 

A self-consistent calculation of the screening current of the NS 
structure was performed taking into account the contribution from the 
Andreev levels. The analysis of the derived expression suggests the 
paramagnetic contribution to current. For example, Fig. 2 illustrates 
the dependence of the current upon the phase (magnetic field). The 
values of the current j to the left of the extremum account for the 
contribution of the Andreev levels. The derivative of this curve with 
respect to the field is proportional to the magnetic susceptibility of 
the NS structure. It is positive ("paramagnetic") in the region of low 
magnetic fields and negative ("diamagnetic") in high fields. 

Similar behavior is observed when the susceptibility of the NS 
structure is measured as a function of temperature in a pre-assigned 
weak magnetic field: it is "paramagnetic" in the region T < and 
"diamagnetic" atT > Tr up to the critical temperature. Temperature 
dependence of magnetic susceptibility in the NS structure at fixed 
magnetic field will be investigated in detail in separate publication. 

In the absence of the proximity effect in the NS structure, when 
the penetrability of the barrier between the S and N metals is small, 
the electrons of the normal metal are reflected specularly from its 
boundaries. In this case the SN structure is a total of two isolated 



22 



subsystems (normal metal and superconductor) placed into a mag- 
netic field. Because of the Meissner effect, diamagnetic current de- 
velops near the superconductor surface. In normal metal, because of 
the Aharonov-Bohm effect, the quantized spectrum of quasiparticles 
is dependent on the magnetic flux through the cross-section of the 
cylinder. The flux generates a paramagnetic contribution to the sus- 
ceptibility whose quasiclassical parameter of the problem l/kpR is 
small. Hence, in the absence of the proximity effect no competition 
is possible between the paramagnetic and diamagnetic contributions 
in the NS structure, and the reentrant effect is unobservable in such 
NS sample. 

To conclude, it should be noted that the explanation proposed in 
this study for the reentrant effect was developed within a model which 
does not assume the electron-electron interaction in the N layer of the 
NS structure. In terms of the free-electron model, a large paramag- 
netic contribution to the susceptibility of the NS structure appears 
in the region of very low temperatures in a weak magnetic field. If we 
increase the thickness d of the pre-assigned normal metal, this would 
lead to a greater number of the Andreev levels no in the potential 
well and affect the solutions of the self-consistent equation for a. As a 
result, the shape of the curve of the paramagnetic susceptibility would 
be slightly "deformed" though its qualitative behavior would remain 
the same. 
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Appendix 

Let us calculate the integral taken over x in Eq. ( ITHl) 

Xo 



I dxx'^Olan — bx — t\J\ + x' , . , , 

J= / , ' :. Al 

y (a, - 6x) V(an - &a:)2 - e2(l + x'^) 

After introducing the notation = a^/^, we can see that the 
function in front of the radical in the denominator has a singularity 
at the point x = a^. Besides, as was noted in the text, the integrand 
has singularities at the points xi, x^- 

Integral ( lAll) can be written as a sum of four integrals 

J = J dx . . . = lirn < j dx . . . + j dx . . . + j dx . . . + j dx . 

V Xi+e aQ+e X2+s 

It is obvious that the presence of the ^-function makes the second 
and the third integrals equal to zero. We first calculate the integral 

x\—e 

Ji = ilim / ^"^"^^ (A2) 



where A = }? — ^ , B = — 2a„6, C = — e^. On substituting the 
variable an — x = l/l the indefinite integral becomes f dtt{ant-i) 

where a = — (1 + Q^o)^^' = 2aoe'^. It can be calculated by the 
method of undetermined coefficients: 



/ 



dtf{t) 



f dt 



if f{i) is the polinomial to power n. Although the calculation is te- 
dious, it is actually simple. The coefficients Ai^ A2, A^ and A4 are 
readily found as: 

2al -3 + 5al + al/2 
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^ ^l{all2 - 1) ^o(2 - all2) 
' e2ao(l + a^)2 (1 + «o)' ' 
It is seen that the coefficients have different orders of e^^-magnitude: 
Ai,A2,A4 ~ e^^, ~ e~4. Finally, we have to calculate six integrals 



t-i—e 



Ji = ilim J dt|2^ity^+^y^ + ^iatV/R(t)+ 

(A3) 



where R{t) = + [3t + A and the designations to = ^5 — £^ = (o^o — 
xi + e)^^ are introduced. All the six indefinite integrals in expression 
31) can be calculated accurately After substituting the limits 



of integration, integrals 1, 2, 3, 4 and 5 are bounded above on energy, 
which is due to the term ^/ R{t) ~ y^a^ — e^/ag, i.e. 9{an — e). Taking 
into account the determined coefficients Ai [i = 1,2,3,4), we can 
obtain the final expression for Ji. 

^ ^ r bjal + 1/3) _ bil + al/6) 2a^6 
'~ 6U3e2(l + a2)2 e2(i + ^2)2 ^ 3e4(i + ^,2)2 + 

6(1 - fqg - a^/6) ^o62(ag/2-l) ao(2-ag/2) 

e2(l + a2)3 + 63(1 + ^2)5/2 + £(1 + ^2)7/2 ^ ^ 

6(a2/2-l) (2-a2/2) -| 

£2(1 + ^2)2 (1 + 0,2)35 J- 

Of all the terms in ( 1A4[ ) , the most significant contribution is made by 
the third term because there is a factor in the numerator of the 
integral over the energy in Eq. (1171) . The contributions of the other 
terms are negligible. We thus obtain the estimate 

2al 1 

^ 36(l + a2)2^- ^^^) 



A similar calculation of the integral 
J4 = — lim 



Xo 

1 f dxx' 



£-0 J {x- ao)WAx^ + Bx + C 

X2+e 



25 



gives a contribution, which is identical in the order of magnitude with 
(IA5I ) . As a result, we obtain the J estimate present in the text, Eq. 

The author is sincerely grateful to A.N.Omelyanchouk for helpful 
discussions and support, to S.I.Shevchenko for valuable comments. 
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